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After summarizing the current theoretical status of the four-body decay B —>■ K*{-^ fi~, 

we apply the formalism of spin-symmetries to the full angular distribution, including the S-wave part 
involving a broad scalar resonance Kq. While we recover in the P-wave sector the known relation 
between the angular observables pI'\ we find in the S-wave sector two new relations connecting 
the coefficients of the S-wave angular distribution and reducing the number of independent S-wave 
observables from six to four. Included in the experimental data analysis, these relations can help to 
reduce the background from S-wave pollution. We further point out the discriminative power of the 
maximum of the angular observable P 2 as a charm-loop Insensitive probe of right-handed currents. 
Moreover, we show that in absence of right-handed currents the angular observables Pi and P 5 fulfill 
the relation Pi = PPI, at the position where P 2 reaches its maximum. 


I. INTRODUCTION 

Rare B decays constitute one of the cornerstones in the search for physics beyond the Standard Model (SM). Among 
them, the semileptonic mode B —> K*{-^ ii~ represents a particularly interesting channel as the measurement 

of the 4-body angular distribution provides a plethora of information which can be used to probe and discriminate 
different scenarios of New Physics (NP). In 2013, LHCb presented results of the measurement of an optimized set 
{Pi'^} of angular observables [THS] based on lfb“^ data. These observables are constructed in such a way that, to 
leading order in the strong coupling constant as and in the large-recoil expansion, non-perturbative form factors cancel 
in the region of low squared invariant mass of the dilepton pair, a unique and powerful feature in the hadronic 
environment. 

Experimental data showed several interesting tensions with respect to SM expectations : Most striking is the 4 a 
anomaly^ encountered in the observable P 5 ID in the bin [4.3,8.68] GeV^. The observable P 2 [HIS] further displayed 
a 2.9cr deviation in the q^-bin [2,4.3] GeV^. The position of its zero {qq = 4.9 ± 0.9 GeV^), which is identical to 
the zero of the forward-backward asymmetry Ape, is in agreement with the SM prediction ~ 4 GeV^ but allows 
for higher values. It is remarkable that all these deviations point to the same negative NP contribution to the 
Wilson coefficient of the semileptonic operator Og, possibly accompanied by a NP contribution to the Wilson 
coefficient of the magnetic operator Or. New Physics contributions to the Wilson coefficient Cio, and, in particular, 
to the coefficients g of the chirality-flipped operators are consistent with zero already at 1 a. The full pattern, 
first pointed out in Ref. [0] and obtained using all available experimental bins in B ^ K *together with data 
on R —)• K*"f, B —>• As 7 , B —)• and Bs —>■ is given by the 1 a ranges 

e [-1.6,-0.9], G [-0.05,-0.01], CfgP G [-0.4,1.0], 

CfP G [-0.2,0.8], CfP G [-0.04,0.02], G [-0.4,0.4], (1) 

where the mild preference for a positive is mainly driven by Bg —> lP~ data. 

The large negative NP contribution to Cg was independently confirmed later on by other groups, using different ob¬ 
servables Si [HKin] (relying on the single large-recoil bin [1,6] GeV^ and low recoil data), different statistical approaches 
Em or form factor input from lattice [12]. Although it had been shown in Refs. mm that a large -I- Cg < 0 was 
preferred in order to explain the P^ anomaly, the possibility of a substantial positive Cg enforcing C|^^ -I- Cg ~ 0 was 
discussed in Refs. [iniEiij. driven mainly by the 1 fb ^ data m on the charged B decay R+ —>■ yS 'in the region 
of low hadronic recoil. The situation has become more coherent recently as the latest 3 fb“^ data on R+ —^ Ar+/i+/r“ 
and R° —?► provided by LHCb [16] is also in good agreement with the solution Cg^ + Cg < 0 (5, both in 

the region of large as well as low hadronic recoil [13 HI]. The three modes thus seem to point to a consistent overall 


^ In Ref. [3 this discrepancy is quoted as a 3.7 cr tension between the experimental result and the 68.3% confidence level of the theoretical 
prediction, while we have quoted the tension between the experimental result and the theoretical central value. Note also that using the 
updated predictions [8] for all observables, including parametric and form factor errors, factorizable power corrections together with an 
estimate of non-factorizable ones and charm-loop effects, the tensions with data, albeit slightly reduced, are still clearly present. 



2 


picture of NP in agreement with the pattern given by Eq. Q. Moreover, under the assumption that NP affects only 
muons but not electrons, also the 2.6 a deviation measured by LHCb [19] in the observable 

_ Br(i3+ ^ K+^i+^l-) 

^ Br(B+ ^ K+e+e-) ^ ’ 

can be explained within the same scenario [2ni - l22| . In order to be able to draw solid conclusions and to see how this 
pattern evolves, it will be crucial to know the 3fb“^ data on the observables in B ^ K*. 

In parallel, the question has been raised if the observed discrepancies between data and SM predictions could be 
attributed to non-perturbative QCD effects [23], even though hadronic form factors enter optimized observables Pj'^ 
only at order ag or through corrections breaking the large-recoil symmetries (factorisable power corrections). There 
exist two different approaches to account for factorisable power corrections: they can either be calculated (under certain 
modelling assumptions) within a non-perturbative framework like light-cone sum-rules (LCSR) [lOl 125] , or they can be 
estimated exclusively on the basis of dimensional arguments and fundamental model-independent relations [U [23l [28| . 
While the first method with full correlations among the form factors is suitable in order to extract the maximal 
information from a particular non-perturbative calculation, the second option in which correlations are included 
via large-recoil symmetry relations reduces the dependence on non-perturbative input to a minimum. The two 
approaches are thus complementary and, because the large-recoil symmetries are expected to be the dominant source 
of correlations, they should give similar results. Indeed, the resulting uncertainties obtained with the first method in 
Ref. (TD] and with the second method in Ref. [5| are of the same order of magnitude (see also Ref. jUj). Both these 
analyses find hadronic uncertainties from form factors to be under control^. 

As a different explanation of the anomaly, the possibility of a large non-perturbative charm-loop contribution has 
been proposed |29j . requiring a huge correction with respect to theory predictions within the factorization approxi¬ 
mation. The discussion in Ref. |29j relies on two model-dependent assumptions: First that the resonance structure 
obtained from a fit to high-g^ data on the scalar mode —>■ K'^fjT can directly be transferred to the vector mode 

B —>■ K*pL~, and second that it can be extrapolated to low values of . The only existing calculation [30] seems to 
be in contradiction with the low-g^ scenarios of B —>■ K *obtained with this ansatz in Ref. [55] as it finds a much 
smaller size for the charm loop and, moreover, the opposite sign for its contribution in B^ —>■ K'^p~ as compared 
to B —>■ K*p^p~^ (contrary to the assumption in Ref. [29]). Furthermore, if the 2.6cr deviation in the observable 
Rk persists, it poses a serious problem for the charm-loop or any other low-energy QCD explanation which cannot 
generate effects violating lepton-flavor universality. Also the observable P 2 in B —> K* p'^ p'^ can be instrumental in 
testing the charm-loop hypothesis proposed in Ref. [55] (see also [55]b 

While the polluting effects from non-perturbative QCD have been studied in detail in the literature, less attention 
has been paid to the so-called S-wave pollution, generated by the background decay B — Kq{-^ K'k)p'^ p~ where Kq 
is a broad scalar resonance. In Ref. m a detailed and complete calculation of the S-wave background was performed 
and it was concluded that any observable will unavoidably suffer from its pollution. While this conclusion is correct 
in the case of uniangular distributions, it does not apply to full or folded distributions where the P- and the S-wave 
parts can be separated according to their different angular dependence. As shown in Ref. |32| . S-wave pollution can 
be avoided for the Pj^'^ observables if folded distributions are used instead of uniangular ones. A discussion of the 
experimental implications of the S-wave contribution was presented in Ref. [33] (see also Ref. |33|b 

Experimental analyses of B —> K* p'^p~ rely on theoretical information regarding the S-wave background. To 
this end, a set of model-independent bounds on the coefficients of the S-wave part of the angular distribution was 
presented in Ref. [3, derived from application of the Cauchy-Schwarz inequality. On the other hand, it was shown 
in Refs. [3131] that the coefficients of the P-wave part are not independent parameters but that they are correlated 
through the spin-symmetry of the angular distribution. In this work we transfer this idea to the S-wave sector. We 
derive two relations which effectively reduce the number of free coefficients of the S-wave distribution from six to four. 
It is expected that the inclusion of these relations into the data analysis can help to further improve the background 


^ Refs. |23l I28| . on the other hand, quote much larger uncertainties. One of the reasons for that has been identified in Ref. the 
decomposition of a form factor into a leading-order part and a OfA/m;,) power correction is not unique but (as in any fixed-order 
calculation) introduces a scheme dependence of observables at neglected higher orders in A/mi,. As the observables are effectively 
calculated at leading order ( 0 (A/m 5 ) effects are not calculated but only estimated), they exhibit a scheme dependence at 0{A/mi,) 
implying a ~ 100% scheme dependence of power corrections. In order to ensure predictivity of the method, it is hence crucial to exploit 
the freedom of choosing a scheme to minimize the impact of the unknown power corrections on the relevant observables (in the same 
way as in a fixed-order calculation in quantum field theory the renormalization scheme is chosen such that neglected higher orders do 
not spoil the perturbative expansion). It was demonstrated in Ref. that the sensitivity to factorizable power corrections of the key 
observables like is significantly reduced if a different scheme is chosen than the one employed in Refs. |23lI28| . 
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estimation. We illustrate the effect of the correlations for the ratio of the S-wave observables Ag and and study 
implications at the position = qf of the maximum of the observable P 2 - Moreover, we point out a relation between 
P 4 and P 5 at q^ = q\ and suggest to use the maximum of P 2 as a golden observable to probe right-handed currents 
(for an explicit model generating right-handed currents see e.g. Ref. [5B]b 

The outline of this paper is the following: In Sec. [h] we discuss the spin-symmetry of the differential B —>■ 
decay rate and determine the number of independent observables in the P- and in the S-wave sector. In Sec. |III| 
we derive the resulting symmetry relations. Their phenomenological consequences are discussed in Sec. HYl First, 
we study the discriminating power of the maximum of P 2 as a test for right-handed currents, then we determine a 
relation between P 4 and Pg at the position of the maximum of P 2 , and finally we investigate constraints on the S-wave 
observables A^'g and derive simple relations among them at the position of the maximum and the zero of P 2 . Sec. ^ 
contains our conclusions. In Appendix A we present an explicit example of how to use the freedom introduced by the 
symmetries to fix a possible convention for the amplitudes, while Appendix B contains details of the derivation of the 
bound on A 5 . 


II. SPIN SYMMETRY OF THE DIFFERENTIAL DECAY RATE 


The differential decay rate of the full four-body decay B —>■ KTr£~^£~ receives contributions from the P-wave decay 
B —>■ Ar*(—>■ Ktt)£~^£~ as well as from the S-wave decay B —>■ A'q(—>■ K'k)£'^£~ with Kq being a broad scalar resonance. 
It can thus be decomposed into a P-wave and an S-wave part, 


dq^ dcos 6 k dcos 9^ dcj) 


Wp + Ws, 


(3) 


with Wp containing the pure P-wave contribution and Ws containing the contributions from pure S-wave exchange 
as well as from S-P interference. Here, q^ denotes the square of the invariant mass of the lepton pair and rriK-K the 
invariant mass of the Ktt system. Further, 0^, 9k are the angles describing the relative directions of flight of the final- 
state particles, while 4> is the angle between the dilepton and the dimeson plane (see Ref. [35] for exact definitions). 
Angular momentum conservation dictates the dependence of Wp and Ws on 9^, 9 k, (j) to be 


and 


Wp = 7 ^ [Jis sin^ 9 k + Jic cos^ 9 k + {J^s sin^ 9 k + J^c cos^ 9 k) cos 29i 

3z7r 

+J 3 sin^ 9k sin^ 9i cos 2(f> + J 4 sin 29k sin 29i cos cj) + J 3 sin 29k sin 9i cos 4> 

+{Jes sin^ 9 k + Jqc cos^ 9 k) cos 9i -\- J 7 sin 29 k sin 9i sin ())-{- J 3 sin 29 k sin 29i sin (j) 
+Jq sin^ 9k sin^ 9i sin 2^] 


= 4^1 


Jia + Jib COS 9k + (TIa + '^26 COS 9k) COS 29 1 -|- J 4 sin 9k sin 29( cos (j) 
+J 3 sin 9k sin 9i cos (j) + Jj sin 9k sin 9e sin (p + J 3 sin 9k sin 29 1 sin p 


(4) 


(5) 


The coefficients Ji and Ji are functions of q^ and 

If not explicitly stated otherwise, we will neglect lepton masses in the following. Then, the decays B —>■ K*£~^£~ 
and B —)■ Arg£+£“ are described by six complex amplitudes A^p^^ and two complex amplitudes Ag^’^, respectively, 
where the upper index L,R refers to the chirality of the outgoing lepton current, while in the case of the P-wave 
the lower index ||,-L,0 indicates the helicity of the Ar*-meson^. These amplitudes are multiplied by a Breit-Wigner 
propagator BWi{m\^) with i = K*,Kq describing the propagation of the K* and Kq meson, respectively. For the 
exact form of the Breit-Wigner functions BWi{m\^) we refer to Ref. m- 

Since the final state is summed over the spins of the leptons, the obervables Ji and Ji are exclusively described in 
terms of spin-summed squared amplitudes of the form This pattern suggests to combine left- and 


^ In the case of non-vanishing lepton masses and of scalar operators coupling to the lepton pair, two additional amplitudes At and As 
have to be included. 

Interferences of different K* and Kq helicities i ^ j contribute, as these particles only appear as unobserved intermediate states. 
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right-handed amplitudes to two-component complex vectors: 


/ Af^BWp 


nj_ = 


A^BWp 

-A^*BW^ 


_ f A^BWp 

~ \A§*BW^ 


ns = 


f A'o^BWs 

\A'g^*BW^ 


( 6 ) 


Using this notation, the observables Ji and Ji can be expressed in terms of scalar products n\nj of these vectors. 
Neglecting lepton masses and presence of scalars we find 


Jis = ^ (|n_Lp-h |n||p) , 

J2c = —\nof, 

J 5 = v^Re(nJn_L), 

Js = -^Im(nJn_L), 


Jlc — I'^ol 


•^3=2 “ 

Jqs = 2Re(n^ U||), 
Jg = -Im(n^U||) 


J2s = ^ (|n_Lp -h \n 

Ji = ;^R-e(4n||) 
J 7 = —'\/ 2 Im(nQ ny 
Jqc — 0 5 




(7) 


and 


rc — _ jc 
^la — ^ 2 a 


rc — _ T<J 
^Ib — ^ 2 b 


= -{\A',^\^ + \A',^nBWs\^ = -\ns\^, 

= ^ySRe [{A'q^AI^* + A'q^A^*)BWsBW^] = ^y3Re(4no), 


3 /3 


Ji = TA/ 7:^-6 


4 V 2 


(Tli 




L* 


A'o^Al^*)BWsBW^ 


- Ao'^A^*)BWsBW*p 
Jr = ^y|lm[«Af-A'«Af)RlUsRfU7 
iVflm [«A 


Js = 


Js = 


+ A'J^Af’)BWsBWp] 


-Re(n^n||), 


-Re(n^ n_L), 


^Im(n| ns), 


^lm(n^ns). 


( 8 ) 


The fact that the Ji and Ji observables involve a sum over the spins of the leptons implies that they are not sensitive 
to the full information contained in the helicity amplitudes Ag^’^. This can be easily seen from the notation 

in terms of the vectors ni. As the Ji and Ji observables are scalar products of the ni, they are invariant under a U(2) 
rotation of these vectors. It is thus impossible to fully reconstruct the amplitudes from the Ji, Ji observables alone. 
If one wishes to extract the A|j'’_[^g, Ag^’^ from experiment, it is mandatory to fix a convention which resolves the 
ambiguity related to the U{2) symmetry. A possible choice is presented in Appendix A. 

The number of independent observables that can be constructed from complex amplitudes is given by 2n^. In 
presence of a symmetry S with Ugen generators, there exist 

r^Ohs ~ 271^ TTgen (9) 

independent observables which respect the symmetry S. The 17(2) spin symmetry of the Ji and Ji observables with 
rigen = 4 generators thus leads to the following consequences: 

• In the P-wave sector there are = 2 • 6 — 4 = 8 independent observables. This observation implies the 
existence of a relation between the 9 non-trivially different P-wave coefficients Ji. The corresponding relation 
has been derived in Ref. [35] and its phenomenological consequences have been discussed in Ref. m- 

• In the S-wave sector there are = 2 • 8 — 4 — = 4 additional observables. This observation implies the 

existence of two additional relations among the 6 S-wave coefficients Ji and the P-wave coefficients Ji. These 
relations will be derived in the following section. 

This parameter counting implies that a basis in the P-wave sector consists of 8 independent observables, like the 
basis {r'j Afb or Fl, Pi, P 2 , J 3 , J 4 , J 5 , Jg} proposed in Ref. [ 3 . In particular, the observables of this basis are not 
related among each other through a symmetry, but they are connected to the observable Pg. In the S-wave sec¬ 
tor, a basis consists of 4 independ ent observables. This means that from the complete set of S-wave observables 
{P 5 , A 5 , Ag, A|, Ag, A|} (see Sec. Ill for their definition) a subset of four has to be chosen as basis, while the 


remaining two are obtained from symmetry relations. 
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III. P-WAVE AND S-WAVE SYMMETRY RELATIONS 


The observables Ji and Ji can be expressed in terms of scalar products n\nj. Since ny and n_L span the space of 
complex 2 -component vectors, the other two vectors can be expressed as linear combinations of the former: 


Ui = ainy -I- 5in_L, i = 0,S. 

Contracting with ny and n_L we get a system of linear equations 

= Oilnyp-I-&i(n|n_L), 
n\ni = ai(n^ny)-I-6i|n_Lp, 

which can easily be solved for aijbf. 


Ui = 


|n_Lp(n|ni) - (n|n_L)(n^ni) 

|ny|2|n_L|2 - |n^ny|2 


nyp(n^ni) - (n^ny)(nj|ni) 

|ny|2|n_Lp - InlnyP 


h = 


( 10 ) 


( 11 ) 


( 12 ) 


Using the decomposition (10 1 of no,ns in terms of ny,n_L to calculate the scalar products |n.op, \ns\'^,nQns, one finds 

= ai{nlnn)+h{nln±), {i = 0,S) 

nlns = as{nlni\)+ bs{nln±). (13) 

Reexpressed in terms of the coefficients Ji, Ji of the angular distribution, this gives the three symmetry relations^: 
J2c [l6j|s — (4j| -I- j|g -I-4Jg)] = 4 [J6 s(J 4>/5 + <^7>^8) + •^9(<1^5>.l7 — 4 J 4 Jg)] 


— 2 [(2 Jgs -I- J 3 ) ( 4 J 4 -1- J 7 ) -I- (2 Jgs — J 3 ) ( J| + 4J|)] , 


(14) 


--JfJl6J|,-(4j| +J|,+4J2)] 


= 4 


-2 


jQs{JiJb + JlJs) + J9{J5J7 — 4 J 4 J 8 ) 

( 2 J 2 S + J 3 ) ^ 4 J 4 -\- J^^ + ( 2 J 2 S — J 3 ) [jf + 4j| 


(15) 


2Jf, [16J|, - (4J| + Jl + 4J2)] = -4 Jos{Jih + + Jik + JbJi) + + J 7 J 5 - 4 J 4 Jg - 4J8J 4 ) 


+4 (2J2s -I- J3)(4J4 J 4 -|- J 7 J 7 ) -I- ( 2 J 2 S — Jo){JbJb + 4J8Jg) 


(16) 


Eq. (14) had already been derived in Ref. 


, determining explicitly the amplitudes in terms of the Ji coefficients 
after fixing a “gauge convention” (see Appendix A for a possible gauge condition). Here, it has been obtained in a 
“gauge-independent” way. As a cross-check, we have also rederived Eqs. (15),(16) following the alternative procedure 
of Ref. [35]. Of the two relations involving S-wave parameters, eq. (15) and eq. (16), the first one is quadratic in the 


Ji while the second one is linear. It is interesting to note that relation (15) for the S-wave coefficients 74 , 5 , 7,8 has the 
same structure as the well-known relation (14) for the P-wave coefficients 74 , 5 , 7 , 8 , and further that the combination 
of the three equations for 72c — |7f^ =F 27^^ has exactly the same structure as Eq.(14) substituting 7^ —>■ 7^ ± Ji for 
i = 4, 5, 7, 8. 

The S-wave observables are defined as 


A 5 = - - 


3r: 




full 7 hfull 

Ji Ji 


CP _ 8 7ff, 7f^ 
- 0— -^7- 


3r 


full 


full 


3ru + r: 


AbCP = 4 


4 Ji Ji 


full 


3r: 


(17) 


full 


full 


5 


The same results are obtained if instead of {n||,np} a different subset {ni,nj} is chosen as basis and the derivation is adjusted 
accordingly. In particular, the stated results are valid also for values of for which nii and n± become aligned. 
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where Ji and denote the corresponding angular coefficients and differential decay width for the CP-conjugated 
decays B —)• and B —)• The total differential decay width Fj^jj is given by 


F' — F' 

^ full — ^ K* 


F' 


KX 


(18) 


where in the limit of massless leptons 


F)f. — 4J2 s — J2c, 


^ if * — 3 ''ia- 


(19) 


Expressing Eqs. (14|-(16) in terms of the S-wave observables and 


Fs = 


^0 _ ^0 

F' -I- F' 

^ full ^ full 


and the P-wave observables and Ft (as defined e.g. in Ref. 0 ) we obtain 

kL [4 - P? - 4/f - 4P|] = -4P2 [PiPi + PliP^] - 4P3 [Ps^e " ^^4^8] 

+ (fcT + Pi) [(P.^)" + (P^)"] + {kT - Pi) [(P^)" + {P^f 
ksFrFsil - Fs) [4 - pf _ 4P| - 4P|] = -^P 2 + ^Pg - AA^As] 


As 


Fj' 


1 - Ft 


[4 _ p2 _ 4P| - 4 P 32 ] = -4P2 [Pi^l + 2P'4 - 2P'A| - P'AJ] 


with 


( 20 ) 


( 21 ) 


+ 3 (^t + Pi) [4(^1)^ + (A's)'^] + -{kr - Pi) [(^s)^ + 4(4l|)^] , (22) 


+ 4 P 3 [P^Al - P;,A% - 2PiA% + 2P'4] 

+2(At + Pi) [2P[A% - P'^Al] + 2(fcr - Pi) [P',A% - 2P'4l|] , (23) 


ks = kr = ks = 1 . 


These relations are valid up to terms which are quadratic in the CP-violating parameters Ag 
F!^^. Exact versions of the equations can be obtained by the replacements 


(dcp 


p(') _^ p(') _ pi') 


3 (/)CP 


fc, = i + e;^^/p, 


A^ 

{i = L,T,S), 


(i) CP 




(i) 


■A 


4 ) CP 


(24) 


FF, i^WCP 


(25) 


or 


p(') _V p{'J _ p{'J _ p 

h ^ k = l-F^^/F, 


{') _ p(') 


3 (/)CP 


.(f)CP Ui) _ -0) 4 b)CP 

' -^S — "^5 ’ 


{i = L,T,S). 


(26) 


In the form the equations are displayed, lepton masses are neglected. For the P-wave observables, the full lepton-mass 
dependence can easily be restored by the replacements P 2 —>■ PP 2 , P 5 —t fiP^ and Pg —/3Pg with /3 = ^1 — 4 to^ / q'^. 
In the following we will typically suppress factors of /3 « 1 and only restore them in final results. For the S-wave 
observables, given their poor experimental precision, we will neglect any terms suppressed by small lepton masses 
throughout the paper. 


Note that Eq. (21) is equivalent to Eq. (4) of Ref. [37], while Eqs. (22),(23) involving the S-wave parameters 
constitute the main result of the present work. The information contained in the two S-wave relations is twofold. On 
the one hand, they can be used to obtain independent bounds on the five S-wave observables As^A‘g^A%, As 

we will show, the resulting bounds are equivalent to the ones derived from the Cauchy-Schwarz inequality in Ref. |S| . 
On the other hand, the equations relate the six S-wave observables Ag,Ag and Fs to each other, reducing 

the number of independent observables effectively from six to four. These correlations should thus be implemented 
in the experimental data analysis in order to improve the background estimation. 
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IV. PHENOMENOLOGICAL IMPLICATIONS 
A. Connecting Pi and P2: The maximum of P2 as a test for the presence of RH currents 


Before discussing the phenomenological consequences of Eqs. (21 )-(23), let us first have a closer look at the observable 

(27) 


x = ki,-P't- 4/?2p| - 4P 


appearing on the left-hand side of these equations. In eq. (27) we have reinstalled the dependence on the lepton 


mass by means of the parameter /3 = — 4 to |/. Expressing kT,Pi, 2,3 in terms of ny and n_L (and the respective 

vectors hy and n± parametrising the CP conjugated amplitude), it can easily be shown that a; > 0 up to terms 
quadratic in the CP-violating observables® From this observation the upper bounds |Pi| < 1, 

IC 2 I < l/(2/3) and jPaj < 1/2 can be read off immediately. On the other hand, it can be concluded that, if one of the 
three observables P’ 1 , 2,3 saturates its bound at a point = qf, the other two observables have to vanish at this point. 
The experimental result (P 2 ) [ 2 . 4 . 3 ] = O.SOj/o]™ indeed suggests a quasi-saturation^ of the bound for the observable P 2 
in the bin [2,4.3] GeV^. Depending on how this result evolves with the new data, the correlation with Pi via the 
positivity condition a; > 0 could be useful to constrain the less precisely measured observable Pi in the respective bin. 

In order to study the information encoded in the maximum of P 2 and the relation with the observable Pi in more 
detail®, let us have a look at the expressions of these observables in terms of the vectors n± and ny: 


Pi = 


KP-hiP 




hiP’ 


P 2 = 


2/3 


1 - 


(n_L - ny)l(n_L - ny) 


+ \n\\ 


(28) 


Obviously, P 2 reaches the extreme value I/(2/I) at the position q^ of its maximum if and only if n±{qf) = ny(g^), i.e. 
if A^{q‘l) = A^{q\) and A^{q\) = —A^{q‘l). At leading order, the second of these two conditions is automatically 
fulfilled in the absence of right-handed currents C/ = Cg = Cjg = 0, while the first condition is fulhlled in this case 
(and neglecting the small ImCI® entering P 2 quadratically) for 


2 2mbMBCf 

9i = 


Cio - ReCfiqf) 


(29) 


From this observation we conclude that any CP-conserving new-physics contribution added to the Wilson-coefficients 
^ 7 , 9.10 will shift the position qf of the maximum of P 2 , while maintaining its height at p^^^ ^ l/(2/3). Compared to 
the SM-prediction qf « 2GeV^, the experimental result (P 2 )[ 2 , 4 . 3 ] = O.SOPq q® prefers a larger value for ql, more to 
the center of the bin [2,4.3] GeV^. This pull to larger g^-values for the position of the maximum of P 2 is consistent 
with the pull to larger g^-values of its zero mentioned in the introduction. From Eq. (29) we see that a larger qf can 
be obtained by a negative NP contribution to Cg, as required by the P 5 anomaly, andTor by a positive contribution 
to CiQ. Notice further that, while it was claimed in Ref. |29] that charm-loop effects might affect the position of the 
zero of P 2 , their impact on the position of the maximum is basically negligible for all scenarios studied in Ref. [29) . 
In general, the maximum of P 2 probes the Wilson coefficient Cl® in a different region in q^ than the P 5 anomaly 
or the zero of P 2 . While a potential NP contribution to C|® is g^-independent and thus induces exactly related 
effects in the three observables, a charm-loop contribution enters C|® as a non-trivial function of q^ which is expected 
to decrease with increasing distance to the cc resonance region. A measurement of the maximum of P 2 can thus 
help to discriminate between NP at high energies and non-perturbative charm effects, and the upcoming data with 
smaller-sized bins will help to determine it more precisely. 

In contrast to Cy.g.ig, a new right-handed contribution to one of the Wilson coefficients C7 g will not only shift 
the position qf of the maximum of P 2 but will also lower its value pushing it below I/(2/3) [5]. At leading order, 

this can be seen from the fact that in the presence of right-handed currents the identity A^ = —A^ does not hold 
anymore for all q'^ so that the two conditions A^^qf) = A^{ql) and A^{ql) = —A^^qf) required for P 2 {qi) = l/(2/3) 
cannot be fulhlled at the same point qf. In general, right-handed currents will cause \ni\_{q\)\ ^ |?^||(di)| and thus 
induce a substantial non-vanishing Pi{q\)^ preventing P 2 (q 2 ) from reaching its absolute maximum \/(2j5). 


® The observables x and x constructed from eq. l |27[ l via the replacements l |25| and ( |26[ l, respectively, fulfill x > 0 and x > 0 exactly. 
^ Note that in practice a complete saturation cannot be accomplished due to the finite bin-size. 

® We will assume real Wilson coefficients for all this discussion. 

















FIG. 1: Left: Comparison of the P 2 -curves (central values) in the SM (green) and in two scenarios of New Physics. The 
scenario NP (red) corresponds to = —1.5, the scenario RHC (blue) corresponds to Cg = 1, C[q = 0.4, Gy = 0.06. Dashed 
lines represent the central value for the integrated bin [2,4.3] GeV^ of the respective curve, while the black cross indicates the 
measured value in this bin. Right: The analogous curves for Pi with the black crosses representing the measured values in the 
respective bins. 


In order to illustrate the discriminating power of the bin [2,4.3] GeV^ of P 2 , we show on the left-hand side of Fig. 
the curve of P 2 (central value) in the neighbourhood of its maximum together with the integrated result for three 
different scenarios: the SM, a new physics scenario NP with = —1.5, and a new physics scenario RHC with the 
right-handed currents Cy = 0.06, Cg = 1, C(q = 0.4. In the scenario NP, the maximum of P 2 is not lowered but its 
position is shifted to a higher g^-value leading to a better agreement of the integrated result with the measured value. 
In the scenario RHC, on the other hand, the height of the maximum is lowered resulting in a stronger deviation of 
the integrated result from the measured value compared to the SM case. The scenario RHC has been chosen in such 
a way that the central values for all low-g^ bins of Pi fall within the experimental Icr regions, as demonstrated in the 
plot on the right-hand side of Fig. It thus constitutes an illustrative example of a setup with new right-handed 
currents to which the maximum of P 2 exhibits a stronger sensitivity than the observable Pi. 


B. Relation between P4 and P5 at the position of maximum and at the zero of P2 


Eq. (21 1 is quadratic in the parameters Pl,P^,PQ,Pg. The requirement of real solutions for these observables 


constrains the allowed ranges of possible values. For example, demanding a real solution for P 4 from Eq. (21) implies 


0 < A(Pi) = -4 x(P^)2 - 4x{P^f - 4 [{kr + Pi)P^ - 2P2P^ - 2 P 3 P']' + Axh^kr - Pi), 


(30) 


with X defined in eq. (27) and fulfilling x > 0. Hence, the first three terms in eq. (30) are negative definite and each 


of them has thus to be smaller in absolute value than the positive fourth term. From this observation we can directly 
read off constraints on jPgj and jPgj, while constraints on IP 4 I and jPgj can for example be obtained by considering 
A(Pg). The total set of constraints is given by 

\Pi\ < VkL{kT-Pi), \P^\ < ^VkL{kT + Pi), iPhl < ^VkL{kT-Pi), |H^| < VkUkr + Pi)- (31) 

As before, these bounds (with the reinstalled /3-dependence for P 5 and Pg) are valid up to quadratic terms in CP- 
violating coefficients, while exact versions can be obtained via the replacement rules (25) and (26). The constraints 
are obtained for x > 0 and thus are valid for any except for the single point where x reaches its minimum value 
X = 0. Continuity of the P/ then implies the bounds to be valid also for x = 0. 

In the limit x —>■ 0 the third term in eq. (30) has to vanish in order to render P 4 real. Proceeding in the same way 


for the other A(P 5 g g) we obtain four relations at = qf with x{qf) = 0 : 

[(fcT + Pl)P^-2P2P'-2P3P']^2 = 
[{kr - Pi)P^ - 2 P 2 P' + 2 P 3 P 1]^2 = 
[(fcT + Rl)Hi-2P2P^ + 2P3P'],2 = 
[(fcT-Pl)P^-2P2P(-2P3P']^2 = 


0 , 

0 , 

0 , 

0 . 


( 32 ) 
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Neglecting PsP^ g, ^ -P 2 ^ 4,5 and including the /3-factor for P^, the last two equations reduce to 


PM) = 




kr — Pi 
kT P^ 


(33) 


91 


This relation is valid at the zero qf of x where P 2 = \/k^ — P^/2I3. For Pi ^ 1, which is an excellent approximation 
in the absence of new right-handed currents, ql coincides with the position of the maximum « kT/{2j3) of P 2 , 

and Eq. (33) becomes 


pM) = MDpMi)- 


(34) 


While Eq. (331 is model-independent, Eq. (34) only applies if there are no new right-handed currents. Its experimental 


validation therefore provides a test on the size of right-handed currents. 

An analogous relation between P 4 and P 5 at the position q^ = q^ of the zero of P 2 was derived and discussed in 
Ref. [37]. We reproduce it here for completeness. Dropping quadratic terms in P 3 , Pg.s and in the P^^ it reads 


[P4'^ + /3^P^\. = 1-77(90^), 


(35) 


where ?7(go) = [Pi + Pi{Pa — P^Ps )]q^ is completely negligible (of order r]{qQ) ~ 10 ^) in the absence of new 
right-handed currents. Let us assume that, as data seem to suggest, the zero Qq of P 2 would be larger than predicted 
by the SM. In this case, Eq. (35) forces the absolute value of ^ 5 ( 90 ) to smaller than in the SM, in agreement with 
the anomaly. 

In Fig. 1^ we show central values of the theory predictions for the two functions P 4 — f3Pg and (Pj)^ -I- l3'^{Pg)‘^ — 1 
for the SM and the new-physics scenario NP with = —1.5. The zeros of the corresponding curves at q^ = q\ and 
q^ = 9oj respectively, demonstrate that the relations (34) and (35) are indeed fulfilled to excellent precision. 


C. Constraints on the and relations at the position of the maximum and the zero of P2 


Eq. (22) is quadratic in the parameters The requirement of real solutions for these observables 

constrains the allowed ranges of possible values. Eollowing the procedure described in Sec. IV B for the P/, we find in 
a completely analogous manner the bounds 


l^sl < \^/^ksFTFs{l-Fs){kT-Pl), 
lA^I < y/iksFTFs{l-Fs){kT-Pi), 


I All < ^?,ksFTFs{l-Fs){kT + Pi), 

I All < ]^^^ksFTFs{l-Fs){kT + Pi). (36) 


Combining further the Eqs. (21)-(23), one obtains a similar bound on A 5 (see Appendix B for a detailed derivation): 

l^sl < 2V3fcLfcsi"s(l-Ps)(l-PT) (37) 
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qHGeY^) 

Illustration of the constraints on A 


«"(GeV0 



2 4 6 

qHGeY^) 


FIG. 3: 

= —1.5 (right two plots). Blue bands correspond to the uncorrelated bounds from eq. (361. Dashed lines illustrate the 


g obtained from relation (221 in the SM (left two plots) and in the presence of 


correlation between A'^ and A^ obtained from relation (221 for the scenario described in the text (orange for the SM, red for 

cr = - 1 . 5 ). 


The constraints (36 1 and (37) are identical to the ones given in eq. (51) of Ref. [5] up to the Breit-Wigner factor 
F = Zj^/XY present in the latter. The results of Ref. [5] were derived using a different method based on the Cauchy- 
Schwartz inequality. The factor T’ is a consequence of the implicit assumption of a narrow S-wave resonance made in 
Ref. [S], and it has to be replaced by its upper limit Fmax = 1 in the general case. This subtlety has little impact on 
the numerical results given in Ref. [5] as the phenomenological analysis there was performed taking F = 0.9 « 1. We 
further note that once again the stated results in Eqs. (36) and (37) are valid up to quadratic terms in CP-violating 


coefficients, with exact versions being obtained via the replacements (25) and (26). 

Proceeding in an analogous way as in the P-wave case in Sec. |IVB[ we find also for the S-wave parameters relations 
at the position = q\ of the zero of the observable x. The corresponding equations read 


[{kr + Pi)Al - AP^Al + 2P:iA%]^, = 0 , 
[{kT-Pi)A%-P2Al-2P^A%]^, = 0 , 
[{kTFPi)A%-P2A%-2PsA%]^, = 0 , 
[{kr - Pi)A% - ^P2A% + 2PsAl]^, = 0, 


and simplify to 


2 A'^‘-~‘^ 


(Qt) = 


A 


5 I ^T-Pl 

kr + Pi 


and 


Aliqi) = 


<11 


2y4^ 


kr — Pi 
kx + Pi 


(38) 


(39) 


<11 


under the assumption of PsAg P^A^g. For Pi ^ 1, one obtains at the position q\ of the maximum of P 2 : 

‘^A%{q\) = A%{q\) and X^g{q\) = 2A%{q\). (40) 


The symmetry relation (22), together with the implicitly contained relations (39),(40) at the zero of x, imposes 
correlations among the A^g implying constraints that go beyond the individual bounds given in Eqs. (36),(37). To illus¬ 
trate this, we assume that a measurement gives ^ Ag^. In this case, the symmetry relation (22) implies a direct 


correlation between Ag and If for example is measured to be ^1^= <xP^ where a = \J2>FtFs{ 1 — Fs)ks/kL, 
is completely fixed to Ag = ^P^. This scenario is illustrated in Fig. for constant P 5 ~ 6 %. The orange dashed 
curves in the plots on the left are obtained for SM values of P 4 5 , while the red dashed curves in the plots on the 
right correspond to the presence of = —1.5 (in addition the SM curve is shown also in the plots on the right to 
visualize the shift between the two curves). If one of the curves is measured for Ag, the corresponding curve for Ag is 
predicted by the symmetry relation, and vice versa. Note that also the blue bands corresponding to the uncorrelated 
bounds from eq. (36) are slightly different in the SM and in the NP case. 


As in the previous section for the P-wave observables, we give also for the S-wave observables simplified versions of 
the symmetry relations at the zero = qQ of P 2 . Neglecting the small P 3 ,Pgg terms, Eqs. (22) and (23) simplify to 


[(4A|2 + A72)(i + P,) + (yl|2 ^ ^ _ pg^pgp^^i _ 


(41) 
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As{ql) 


2VT^{2AUl + Pi)Pi + ^ 1(1 - Pi)Pi) 
V^(l - P?) 


(42) 


V. CONCLUSIONS 

In this article we have exploited the spin symmetry of the angular distribution of the decay B —>■ K*, both 
in the P-wave as well as in the S-wave sector. We have shown that the symmetry reduces the number of independent 
S-wave observables from six to four, implying two non-trivial relations among the observables Fg, Ag, Ag, Ag, Ag 

and which we derived explicitly. The relations allowed us to obtain individual bounds on the A^g'^ which agree 
with the ones determined in Ref. [S] via the Cauchy-Schwartz inequality. However, the constraining power of the 
symmetry relations goes beyond these individual bounds as they correlate the S-wave observables among each other. 
The implementation of these correlations into the experimental data analysis is expected to reduce the background 
from S-wave pollution. As an example, we have shown how for Ag® ^ Ag® the correlations fix Ag from a measurement 
of Ag and Fg (or Ag from a measurement of Ag and Fg) in the whole range of the of the squared dilepton invariant 
mass q^. We further showed that Ag/Ag and Ag/Ag are completely fixed at a point = qf where qf coincides with 
the position of the maximum of the P-wave observable P 2 in the absence of new right-handed currents. 

We also pointed out the strong potential of the maximum of P 2 for probing NP beyong the SM, in particular the 
presence of new right-handed currents, in a region far away from charm resonances. We have shown that a shift of 
the position of the maximum of P 2 compared to its SM expectation, with the height of the maximum kept 

at the SM value l/(2/3), would be a signal of a NP contribution to the SM-like Wilson coefficients Cy, Cg, Cig. A 
maximum value P^^^ < l/(2/3), on the other hand, would detect the presence of new right-handed currents and thus 
complement information from the (currently not very precisely measured) observable Pi. We have further proven 
and illustrated that for Cj = Cg = C(g = 0 the angular observables P 4 and P 5 fulfill PJ = /IPg at the position of 
the maximum of P 2 , so that any deviation from this relation would equally signal the presence of new right-handed 
currents. 
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supported by FPA2011-25948 and the grant 2014 SGR 1450, and in part by the Centro de Excelencia Severo Ochoa 
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Appendix A: Gauge conditions for the amplitudes 


All the angular observables studied by the LHCb experiment are invariant under a U{2) rotation of the vectors 
defined in Eq. ([^. As a consequence, the amplitudes A^'^ cannot be determined unambigously from experiment. In 
order to arrive at a one-to-one correspondence between the experimental observables and the theoretical amplitudes, 
one has to fix a convention which picks for every class of P( 2 )-related amplitudes a certain representant (similar to 
’’fixing the gauge”). One convenient choice that has been proposed and is used by the Imperial group of the LHCb 
experiment [ 88 ] consists in requiring 

ReA^ = 0, ImA^ = 0, ImAg = 0, ImA^ = 0. 


This choice is not unique, several combinations are possible (see Ref. |85j for a different choice). Starting from an 
arbitrary amplitude, one arrives at the above configuration by means of the U ( 2 ) transformation 


■ 0 


cos 0 — sin 0 


coshj 0 — sinhi 0 

0 


sin 0 cos 0 


— sinhi 0 coshi 0 


with 


tan 20 = 2 


ImA^ReAg + {L -n- R) 


ReA^ -I- ImAg tan0 
—ReAg -I- ImA^tanfl ’ 
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IixiAq + tan 9 — (ReA^ — ReAg tan 9) tan 9 
—Re^Q + ReA^ tan 9 + (Im^^ + IuiAq tan 9) tan 9 ’ 

ImA^ + ImA^ tan 9 — (ReA^ — ReA^ tan 9) tan 9 
—ReA;5^ + Re^;^ tan 9 + + ImA^ tan 9) tan 9 


Appendix B: Derivation of the bound on As 


In this appendix we pre sent an e xplicit der ivatio n of the con straint on the S-wave observable Ag given in Eq. (371. 
Combining the relations (21|-(23) as Q^(pT|)+3b^(p^+at>(|2^ with arbitrary real coefficients a,b, one obtains an 
equation for linear combinations aP' ± (2)bA% of P- and S-wave observables which has the same structure as the 
individual P'- and Ag-relations (21) and (|2^: 


Y{a, b) [kl - Pi - 4P| - 4P|] = -4P2 [{aPi + 2bA%){aP^ + bA%) + (aP' - bAl){aP^ - 2bA%)] 

-AP, [{aPi + bAl)iaP' - bAl) - (aPi + 26A|)(P' - 26A|)] 

+ (^T + Pi) [(flPi + 2bA's'f‘ P {clPq — 

+ {kT — Pi) [(aPg + bAg)'^ + {aP^ — 2A|)^] , (43) 


with 


Y[a,b) = a kg + 3b ksFTFg{l — Fg) + abAg 

2 


Fs 


= kg 


b. Ft 


l-Fs 

Prp 


4ikL 1 ~ Ft 


[l2A:ifcs'Ps(l - Fg){l - Ft) - A|] 


(44) 


Requiring A{aP)^ + 2bAg) > 0 in analogy to Eq. (30) in order to ensure that the observable aP^ + 2bAg is real, one 
finds that Y(a, b) > 0. This condition has to be fulfilled for any possible linear combination, i.e. for any value of a, 6, 
which according to Eq. (|44|) enforces Ag to respect the constraint 


|As| < 2^3kgkgFg{l-Fg){l-FT). 


(45) 
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